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I When neutrinos propagate in a background, their gravitational couplings are modified by their weak inter- 

actions with the particles in the background. In a medium that contains electrons but no muons or taons, the 
matter-induced gravitational couplings of neutrinos are different for the various neutrino fiavors, and they must 
be taken into account in describing the phenomena associated with the neutrino oscillations in the presence of 
strong gravitational fields. Here we incorporate those couplings in that description, including also the effects of 
' a magnetic field, and consider the implications that they have for the emission of high energy neutrinos in the 

C^, vicinity of Active Galactic Nuclei. 

Oh: 

r-| ! 1 Introduction 

It is now well known that the interactions of a neutrino with the background particles can influence the neutrino 
■ properties in essential ways as it propagates through a medium. Those effects can have important consequences for 

I various physical phenomena, such as the Mikheyev-Smirnov-Wolfenstein (MSW) mechanism and its variations in the 

context of the solar neutrino problem, the explanation of the large velocities of pulsars in terms of the asymmetric 
emission of neutrinos from the proto neutron star and the propagation of high energy {E > 10^ GeV) neutrinos in 
the vicinity of cores of Active Galactic Nuclei (AGN) §. 

In all these applications, a common ingredient is the fact that the various neutrino flavors interact with the 
background particles differently. While the standard neutrino interactions in the fundamental Lagrangian have a 
universal character, the universality is broken when the effects of the medium are taken into account. 

In the presence of a medium, these breakdown of the universality of the neutrino interactions includes also the 
gravitational ones. In Ref. |^ the effective gravitational vertex of neutrinos was determined by calculating the 
one-loop contribution to the neutrino stress-energy tensor. In the presence of a static gravitational potential, such 
matter-induced interactions lead to an additional contribution to the neutrino dispersion relations, or equivalently 
their indices of refraction, which depend on the gravitational potential and are not the same for all the neutrino 
flavors. As was stressed in that work, those calculations were based on the standard model of particle interactions 
and the commonly accepted linearized theory of gravity with a universal tree-level coupling. The non-universal 
character of the induced gravitational couplings is a consequence of the flavor asymmetry of the background and not 
of an assumed breakdown of universality of the fundamental Lagrangian. 

On the other hand, it has been observed by Pi'riz, Roy and Wudka that high energy neutrinos originating 
from AGN can have resonant spin-flavor transitions due to the combined effects of the gravitational interactions 



1 



and the presence of a large magnetic field. AGN can be a source of high energy neutrinos and, therefore, a good 
understanding of such transitions is useful to estimate the expected fluxes from these objects in the forthcoming 
searches for high energy neutrinos in neutrino telescopes. 

In the calculations of Ref. the matter, magnetic and gravitational effects were considered. However, the latter 
were introduced using the standard tree-level couplings of neutrinos to gravity and, in particular, the effects of the 
medium on the effective gravitational interactions were not taken into account. This amounts to neglect the additional 
matter-induced gravitational contribution to the neutrino dispersion relation. As we have already mentioned, they 
depend on the gravitational potential and, what is more important for the issue of resonant transitions, are not flavor 
symmetric. In contrast, the tree-level gravitational couplings are the same for all the neutrino flavors, including 
the right handed (singlet) neutrinos, and hence should have no effect whatsoever on the phenomenon of neutrino 
oscillations. 

Motivated by these considerations, in this work we take another look at the subject of neutrino oscillations in 
the presence of a gravitational field. Our objective is to pay particular attention to the effects of the matter-induced 
gravitational couplings that we have already mentioned, with a view to their implications in the context of high 
energy neutrino emission from AGN. 

We start in Section || with a brief overview of the main results obtained in Ref. ||] that are needed in subsequent 
sections here. This includes a summary of the formulas obtained there for the dispersion relations for standard 
left-handed neutrinos in the presence of a gravitational potential, taking into account the matter-induced terms in 
the gravitational vertex of the neutrino. In Section ^ we set up in general terms the equations that are relevant for 
treating the phenomenon of resonant neutrino (spin)-flavor transitions in the presence of a gravitational potential. 
The treatment includes the possibility that a magnetic field may also be present, under the assumption that the 
neutrinos have an intrinsic magnetic moment. Using the results of Section |^ as the starting point, in Section ^ we 
study the effect of the resonant transitions on the the determination of the flux of high energy neutrinos emitted by 
AGN, and finally the conclusions are given in Section |l 



2 Matter effects on the neutrino gravitational interactions 

In Ref. we obtained the dispersion relation that is obeyed by a standard left-handed neutrino in the presence of a 
static gravitational potential. In this section we review briefly those results and consider the corresponding ones for 
a right-handed (singlet) neutrino . 

As already indicated, those results were derived in the linearized theory of gravity, in which the metric tensor is 
written as 

gxp = ?7Ap + 2k/iap, (2.1) 

where rjxp is the flat space metric. We then expand the Lagrangian in the presence of gravity in powers of k and 
keep only the first order terms. In this formulation, h\p appears as the graviton field, which is a spin-2 quantum 
field coupled to the stress-energy tensor, whose interactions can be studied in the flat Minkowskian background. The 
Einstcin-Hilbcrt action for pure gravity is given by 

where R is the Ricci scalar, g is the determinant of the matrix gxp, and G is the Newton's constant. Using Eq. (pj]), 
we can verify that this gives the correct kinetic terms for the spin-2 field if we make the identification 



K = V8^. (2.3) 

Then starting from the Dirac Lagrangian for a given fermion / in the presence of gravity, the coupling of the graviton 
field hxp with the fermion field can be written as 

4//)=-K/.^P(x)fi;)(x), (2.4) 
where the stress-energy tensor operator T^y for the fermion field is given by 



K {x) = \ -^(x) [^xdp + ipdx] V(x) + H.c. - Tixpa"{x) . (2.5) 
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Here £q (x) is the Lagrangian for the free Dirac field, which we write in the expUcitly Hermitian form 



if) 



^^j^df^yj + H.c. 



(2.6) 



From Eqs. (2^) and (2^) it follows that the term associated with the gravitational fermion vertex in a Feynman 



diagram is —iKVj^p\ where 



VxPiP^P') = J bx{p+p')p + lp{p + P )x] - ^Vxp [{^ - mj) + if - mf)] 



(2.7) 



It is easy to deduce by inspection that for a left (j/l) or a right {vn) handed massless neutrino, the corresponding 
quantity is given by 
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(2., 



where x = L,R= ^(iTTs)- 

When the effects of the background medium are taken into account, the result is that the matrix element of the 
total stress-energy tensor T^i/(x) between neutrino states, with incoming and outgoing momenta k and k' respectively, 
is given by 

WLM{k')\%AO)WLAk))=ULMk') (l^j;;"'«'(fc,fc') + A|,';^-«)(fc,A;')) UL,i?(fc) . (2.9) 

The quantity A{^t^'"''(fc, k') represents the matter-induced contribution and its calculation to one-loop was the subject 
of Ref. That calculation was performed for the left-handed massless neutrinos, adopting the standard electro- 
weak couplings for the neutrinos and the other particles. The corresponding result for the right-handed (singlet) 
neutrinos, since they have no standard couplings to matter, is 



A(''«) = 



(2.10) 



The particular formulas for A^'i?'^ are not relevant for us here. More important is the result that they impl y fo r 
the dispersion relation that a neutrino obeys in the presence of a static gravitational potential 0^^^* . Using Eq. ( |2.9| ) 
as the starting point, it was shown in Ref. j^] that the dispersion relation of a standard left-handed neutrino that 
propagates with momentum K is given by 



(2.11) 



while for its antiparticle (a right-handed antineutrino) the sign in front of fej^^^t ^'^^ is the opposite. The coefficient 
^mat i^ usual Wolfenstein term while bQ represents the matter-induced gravitational contribution. Following the 
same arguments, it follows that the dispersion relation o;^'' for a sterile neutrino (either left- or right-handed) is 
simply 



2K(t)'' 



(2.12) 



It is therefore clear that the tree-level gravitational neutrino couplings, which are represented in Eqs. (2.11) and 



(2.12) by the 2K(j)'^^* term, do not play any role whatsoever in the phenomenon of neutrino oscillations, since they 
appear as a common factor in the dispersion relations of all the neutrino states. In contrast, the matter- induced 
terms, which are zero for the right-handed singlet neutrinos and do not have the same value for the three left- 
handed neutrino flavors, affect the neutrino oscillations via the MSW mechanismHl. In those situations in which 



the gravitational contribution can be neglected, Eq. (2.11) reduces simply to the Wolfenstein formula. However, in 
environments in which a relatively strong gravitational potential is present, the matter- induced gravitational term is 
important and must be taken into account in any discussion of MSW- type effects under such conditions. This is the 
subject we take up next. 
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3 Oscillations in a gravitational and a magnetic field 



We consider two neutrino families, the first of which we take to be Vei and denote the second one by which can 



be either one of f^^T- or a sterile neutrino Vs- Each family consists of a left-handed neutrino field VaL{i ■ 



and a 



corresponding right-handed partner NaR that we take to be a weak singlet. In addition to the mixing via the mass 
matrix, we allow the possibility that the neutrinos have an intrinsic magnetic moment coupling. In situations in 
which a magnetic field is present, those couplings can induce the (resonant) spin transitions between the two neutrino 
helicity states. We consider two situations separately, according to the type of the magnetic moment coupling that 
they may have. 



Case I Here we assume that the neutrinos have a magnetic moment interaction defined by the term 

L' = -\^tJiabNaR(Tf,uVbLF'"' + H.C. , 



(3.1) 



a,b 



in the Lagrangian. As shown in the Appendix, if the neutrino flavor amplitudes are assembled in the vector 



(3.2) 



then their evolution is governed by the Hamiltonian matrix 



K + h 



mat 



(m+AtB||)^ 



H 



2K 




(3.3) 



2K 



where each entry in this equation is itself a 2 x 2 matrix. In particular, m is the neutrino mass mixing matrix and /x 
is the magnetic moment matrix defined in Eq. (3.1), while K denotes the magnitude of the neutrino momentum and 
Bt are the components of the magnetic field parallel and transverse to respectively. The terms b-^g^^ and be 
are diagonal matrices whose values depend on the composition of the background medium. At this point it is useful 
to introduce the particle and antiparticle momentum distribution functions 



ffjiPf) 



1 



(3.4) 



where the upper and the lower signs hold for the particle and the antiparticle respectively, ^/ is the chemical potential 
of the fermion / and 



The corresponding total number densities are given by 



(2^ 



// 



nj = 2 



d^P 



(2^) 
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(3.5) 



(3.6) 



In terms of these quantities, 



^mat 



be + Y,f ^fbf for 1^, 

EfXfbf 

for Vs, 



for l^^, 



(3.7) 
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and 







for Vf,, Vr 

for i^c, 



(3.8) 



where 



^1 
Jf 



nj), 



V2GF{nf 
-3(n/ - rij) + 



(fiP dF- 



f 



{2nY2E} dEf 



with 



(3.9) 



(3.10) 



Fj ^ A{2Ej - mj)ifj - fj) . 

In Eqs. (3.7) and (3.S), the sum over / must be made with respect to all the particle species that compose the 
background, and Xf stands for their neutral current couplings. So, for example, for a background composed of 
electrons, neutrons, protons and their antiparticles (/ = e,n,p). 



-Xe = Xp = i - 2 sin^ 



(3.11) 



The explicit formulas for Jf were given in Ref. § for various limiting cases. 

Case II In this case we assume that the magnetic moment interaction term is of the formQ 

1 



L' 



^^^R'^t^'^'^eLF^"' + H.c. 



but that the mass terms in the Lagrangian are diagonal, i.e., 



(3.12) 



(3.13) 



so that there is no mixing in the mass matrix. Denoting by a^,^ and /3p^ the VeL and I'^j^ components of the wave 
function, and writing them in the form 

the Hamiltonian matrix that determines their evolution is given in this case by|^ 



X ■ 



(3.14) 



H = 



' ^ + "mat + "g + 2/f + 2K 



K-b[ 



mat 



2 A' 



2K 



) 



(3.15) 



where b'^^^j^ and 6^°'' (a 

Notice that, in the MSW mechanism involving only standard left-handed neutrinos, the matter effects due to 
the neutral current interactions are the same for all the neutrino species and therefore are not relevant as far as 
the oscillation mechanism is concerned. In contrast, in the situation we have just considered, those neutral current 
contributions have opposite sign for Vf,L and v'^ and are zero for Vg, and therefore must be taken into account. 
In the physical settings that we are considering, they can be important since they involve the nucleons and their 
couplings to gravity. 

Eqs. dJ) and ( Islsl ) set the framework for our consideration of the possible effects that the matter-induced 
gravitational interactions of the neutrinos may have for the phenomenon of (resonant) transitions in the presence of 
a strong gravitational field. As a concrete example, we consider those transitions in the environment of AGN. 



are given in Eqs. ( p. 71 ) and (^ 
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4 Resonant transitions in the vicinity of AGN 



For concreteness, we focus our attention on possible transitions between the and Vt flavors for the following reason. 
Some galaxies have relatively bright centers (as compared to the total photon luminosity of the whole galaxy). High 
energy photons reaching tens of thousands of GeV have been observed from these centers. If Fermi mechanisms are 
responsible for accelerating the electrons in these systems, then protons are also expected to be accelerated by similar 
ones. If this is true, neutrinos are expected to be produced in or/and pp collisions, providing the signature for 
proton acceleration. Tau neutrinos are also to be produced in the same collisions but at highly suppressed levels, 
because both the branching ratios and the production rates for unstable hadrons that decay into tau neutrinos are 
much lower than for muon or electron neutrinos. These two (suppression) factors enter into the calculation of relevant 
intrinsic neutrino fluxes multiplicatively, implying that the intrinsic Ur (and £>,-) hux is relatively rather small. This 
Vt- (and Vr) flux may possibly be somewhat enhanced due to (resonant) transitions between relevant neutrino states. 

Briefly, in p^ collisions, the protons and photons may give rise to high energy Vr (and v-r) mainly through 
p + 7 — > D'g + A° + in addition to producing and mainly through p + 7 ^ A"*" n-\- tt"*" . The production 
cross-section for is essentially up to three orders of magnitude lower than that of A+ production for the relevant 
center of mass energy scale. Moreover the branching ratio of to decay eventually into Vr {^t) is approximately 
two orders of magnitude lower than for A+ to subsequently decay into v^. and through 7r+ . These two suppression 
factors along with the relevant kinematic limits give approximately the ratio of intrinsic fluxes of tau neutrinos and 
electron neutrinos as F^{vr + Vr)/F^{Vf. + Vf.) ~ 10~^. 

In pp collisions, the Vt flux may be obtained through p + p Dg + X. The relatively small cross-section for Dg 
production together with the low branching ratio into Vt implies that the Vr flux in pp collisions is also suppressed 
up to 4 — 5 orders of magnitude relative to Ve and/or fluxes. 

Thus, in both type of collisions, the estimated intrinsic Vr flux from the AGNs is rather small relative to Ve 
and/or fluxes, essentially for the entire relevant neutrino energy range (2 • < A"/GeV < 2 • 10^) [|. However, 
as we point out later in this Section, due to transitions between relevant neutrino flavors, it is possible that F'^(i'r -\- 
Pt) / (Vf, + Vg) 3> 10~^, thus raising the possibility of detection of high energy v.^ in new km^ surface area neutrino 
telescopes within the relevant neutrino energy range 

We flrst briefly describe the matter density and magnetic fleld profiles in the AGN. According to the matter 
density profile in the AGN is given as a function of the distance r from the center by 

PAGN{r) = Pof{x), (4.1) 

where po ^ 1-4 • 10^^^ g cm"'^ (we take the typical luminosity of the AGN to be ~ 10*^ erg s~^) and f{x) = 

a;-2-5(i_o.la;"-3i)-i fora; - (10-100). Here x = r/i?s, where i?s is the Schwarzchild radius i?s =i 3-10" (^l^)ni, 

and in terms of it the gravitational potential is expressed as 0^^^* = — l/2a;. We consider the following magnetic fleld 
profile in the AGN § 

Bagn = Bog{x), (4.2) 

where Bq ^ 1.7 • 10^ G and g{x) = x^^-'^^{\ — O.la;'^ "^^)"'''^. We will use these matter density and magnetic fleld 
profiles in our estimates as an example. 

In what follows, the discussion of resonant transitions in the AGN is divided in two parts, according to whether 
or not the magnetic field terms are important. We discuss the two necessary conditions for resonant transitions 
to occur in these two situations, namely, the level crossing and the adiabaticity at the level crossing, and we also 
comment briefly on the possibility of transitions between active and sterile neutrinos. 



4.1 No magnetic field interactions 



In this situation, the mixing occurs via the mass matrix only. After subtracting a term proportional to the identity 
matrix in Eq. (3.3), the relevant 2x2 Hamiltonian matrix in the (i^ei, Vtl) basis is given by 



H, 



with 



A7?F (Am2^/2X) sin 26*6 

(Am2^/2if) sin 26*6^ -l^Hp 



AHf 



AVf 



2K 



■ cos 29e 



(4.3) 



(4.4) 
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where 9eT is the vacuum mixing angle, Am^^ = — is the mass squared difference between the two neutrino 
mass eigenstates fiL(i = 1,2), and 

AT^i.- V2GfK + Je^/-'"^*). (4.5) 



In contrast to the term proportional to the density rig in Eq. (4.5), the term proportional to the gravitational 
potential depends on the nature of the electron gas. For the case of a classical non-relativistic {(3~^ = T <C me) 
electron background 

Je = -PmeUe. (4.6) 

In the following estimates we use T ~ (1 — 10) eV as an example, which is the value that corresponds to the black 
body (thermal) spectrum temperature for the photons mentioned above This value of T may be considered as 
the average temperature for the remaining background particles as well because of thermal equilibrium (it is mainly 
this ultraviolet bump in photon spectrum that interacts with the Fermi accelerated protons producing the relevant 
unstable hadrons). 



Using the matter density profile given by Eq. (4.1) and for K ^ 7- 10^ GeV, the level crossing condition {AHp = 0) 
can be satisfied for Am^^ ~ (lO"^*' — 10^^^) eV^. Without taking into account the (matter induced) gravity effects, 
the level crossing condition can also be satisfied but for values of Aml^ that are approximately two orders of 
magnitude smaller than the one just quoted. However, the other essential condition, namely, the adiabaticity at the 
level crossing |l0| 

n.^(^y'-'^>i, (4.7) 

\ 2K J TT \AVf\ ~ 

(where AVp = d(AVp)/dr) is not satisfied; that is, Kp ^ 1. Therefore, the matter induced gravity effects do not 
lead to resonant flavor transitions between i/^ and 

4.2 With magnetic field interactions 

There are several possible cases that can be considered, depending on the mass mixing matrix and the magnetic 
moment couplings. We consider the specific situation in which the VeL and i/^j^ fields are mixed by a transition 
magnetic moment interaction without mixing in the mass matrix, as described in Case II in Section |^. In this case, 
after subtracting the term proportional to the identity in Eq. ( p.l5| ) the relevant Hamiltonian matrix in the {I'eL, i^tr) 
basis can be written as 

AHsp -2fiBT 

HsF = U h (4-8) 

-2^Bt -AHsp 



where 
with 

AVsp ~ V2Gp 



A 777^ 

AHsp = AVsp - — ^, (4.9) 

ZK 



(1 + 2X,){n, + Je0^^*) + 2j2Xf{rif + JfC^'^') 

f=n,p 



(4.10) 



— m^g. For simplicity of the notation, in the reminder of this section we take a purely transverse 
magnetic field, so that Bt = B. In the more general case, the formulas given below hold with the replacement 

^For the channel Ufj, — > Ur the matter cfTccts are absent, but there is an interesting possibiUty of vacuum flavor oscillations as supported 
by the recent superkamiokande data . The vacuum flavor oscillation probability expression is given by the familiar expression 

PCl^M ^ I'r) = sin^ 2e^r sin^ I -^^L 

where is the relevant vacuum mixing angle with the corresponding mass squared difference Am^^ . We take Am^^ and sin^ 20fiT 
values (Am?^ ~ 10""^ eV^ and sin^ 26fiT ~ 1) as suggested by recent superkamiokande data concerning the deficit of atmospheric muon 
neutrino Taking the typical distance between the AGN and our galaxy as L ~ 100 Mpc (where 1 pc ~ 3 • 10^® m), it follows that 

P(i/^ i/7y~ 1/2 for 2 ■ 10^ < K/GeV < 2 ■ lO'', thus yielding F{ur + Ur) ~ F(u^ + P^). Similarly, vacuum flavor oscillations between 
and may also take place for certain range of neutrino mixing parameters resulting in P{i's — > i^t) ~ sin^ 29eT-. 
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Figure 1: F{vt- + 

whereas the dotted curve is obtained using Eq. (4.14). In both cases, we use /i ^ 10^ ^^/is 



[cm"2s-isj.-iQ(;V"i] vs K{GcV). The sohd curve is obtained for P - 0.8 [using Eq. (4.12)] 



B ^ Bt- Using then Eq. (^ 



the spin-flavor precession probabihty for constant B and AVgi? is given by 



{2^iBf 



(2^^5)2 + {^HsF? 



Ar 



(4.11) 



Here Ar is the width of the region where B is appreciable. Let us study the various interesting situations that may 
arise from this expression for P{i'e ^ ^^r) under different physical conditions. 



Let us consider first the case in which the contribution of the AVsf term in Eq. (4.11) is small. If Am'^^/2K <C 
2/ii?, the above expression for P^Ve ~^ i^r) then reduces to 



(4.12) 



which is independent of K . Using the Bagn profile given by Eq. (42) and for ^ lO^^^/^s (/is is Bohr magneton) 
10, we obtain P{ve i>r) > 1/2 for Am^^ < 10"^ eV^. For this value of Am^^ (~ 10"^ eV^), the resonant spin- 
flavor precession due to matter induced gravity effects do not take place since the relevant level crossing condition 
is not satisfied. In contrast to the vacuum flavor oscillations mentioned earlier, the value of P in this case can be 
different from 1/2, thus raising the possibility of distinguishing between the two neutrino oscillation mechanisms. 
For Aml^/2K ~ 2/LtS or Aml^/2K > 2fj.B, P{ve Vr) is K dependent and is such that P{ve —> Dr) < 1/2. 

Taking into account now the effects of AVsf, the relevant level crossing condition [AHsf = 0, see Eq. (^^)] can 



be satisfied for Ar 



(10 * — 10 ) eV . This relatively small value of Araj,^ is also interesting in the context of 



supernova explosions ll3| and the Sun[|l^. Furthermore, this value of Amj^ is about 2 orders of magnitude larger 
than the one required for flavor level crossing [see Eq. (4.4)]. This is due to the fact that the nucleon contribution in 
AVsF does not vanish (in contrast to that in AVf), as we emphasized earlier. 
The adiabaticity condition at the level crossing in this case is [fi5[ 



KSF 



2{2iiBf 
\^Vsf\ 



> 1. 



(4.13) 



Using the Bagn proflle given by Eq. (4.2), this condition can be satisfled for ij. ~ 10~ /is, which is an order of 
magnitude smaller than the value that is needed in the context of the Sunp4|. The general expression for P(t'e — >■ i^r), 
including the possible non-adiabatic effects (ksf < 1) is given by 



P{Ve 



0.5- (0.5-Plz)cos26Is, 



(4.14) 
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where tan29B = {—2piB) / AHsf is being evaluated at the neutrino production site and P^z = exp(— ^ksf). The 
expected flux spectrum due to (resonant) spin-flavor precession is calculated as 

F{iyr + Vr) - P{Ve ^ Vr)F''{v, + V,) + [1 - P{v, Dr)^ i>^r + ^r) ■ (4.15) 

A similar expression for F{vr + Vr) due to flavor/spin-flavor oscillations in other channels can be straightforwardly 
obtained with appropriate changes. For i^°(i'e + v^) we use the results from Ref. @ for if > 10^ GeV as an example 
and, since F'^{vr + Vt) / F'^{ve^^ + S'e./j) ^ 1, the explicit form of F'^{i^t + i^t) is not important here. For illustrative 
purposes, in Fig. ^ we display F{h'r + Vt) as a function of neutrino energy for various neutrino spin-flavor transition 
mechanisms, using Eq. ( [4.15[ ). In the resonant spin-flavor precession case, since the adiabaticity condition [Eq. ( 4.13 )] 



is satisfied, we have Plz ~ 0. Also, since the corresponding value of cos26'_b is quite small, P{i'e — > Vt) ~ 1/2 [see 



Eq. (4.14)]. We have verified that a similar behavior of P{ve — > Vt) results if the matter-induced gravity effects are 
omitted, but for a value of Am^ which, for the reason we mentioned earlier, is approximately 2 orders of magnitude 
smaller. 

Another possibility, which we have not considered in detail, is that the transitions occur between active and 
sterile neutrinos. As in the case just discussed, the nucleon contribution in the Hamiltonian is important in this 
case also. In the presence of a magnetic field, the resonant spin-flavor transition may lead to energy dependent 
disappearance/appearance of the active neutrino fluxes even if one uses the stringent astrophysical upper bounds on 
the relevant transition magnetic moment p2[ . 



5 Conclusions 

When neutrinos propagate in a medium composed of electrons and nucleons, their gravitational couplings are modified 
due to the weak interactions with the particles in the background. While the tree level gravitational interactions 
of the neutrinos have no effect on the oscillation phenomena, the matter-induced gravitational couplings depend on 
the background particle number densities and, in the presence of a gravitational potential, they lead to additional 
contributions to the neutrino indices of refraction which are not the same for the various neutrino flavors. 

We have considered the effects that such matter-induced couplings may have on the phenomena associated with 
resonant neutrino oscillations, under the combined presence of a gravitational and a magnetic field. As an example 
of a setting where these effects are relevant, we have studied their infiuence on the determination of the flux of high 
energy neutrinos from the core of Active Galactic Nuclei. In that context we have pointed out that the resonant 
neutrino spin-flavor transitions may take place for Ato^ ~ (10~* — 10~^) eV^ and /i ~ 10~^^/j,b, and due to the 
matter-induced gravity effects the estimated high energy tau neutrino flux is somewhat higher than it would be if 
those effects are not included. We have relied on the perturbative calculation of the matter-induced gravitational 
couplings carried out in Ref. |^ . According to the current interpretation of the observed photon flux from cores of 
AGN's, the neutrino production should take place for values of x around 10 — 100 as we have assumed in the text. In 
this case, the use of those perturbative results is justified and as a consequence the conditions for resonant transitions 
do not depend on the choice of any metric parameters unlike in Ref. ^ . 

We have set aside some possible incoherent effects which could be present for the high energy neutrinos that we 
are considering. The considerations of the present work indicate that more detailed calculations taking them into 
account are worth pursuing, and that they have useful applications in the context of neutrino emission in the vicinity 
of the core of AGN, with interesting implications for km^ surface area high energy neutrino telescopes. 
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Appendix 

A Equation for the amplitudes 



Here we derive Eq. (3.3). We begin by considering one neutrino, with both a left- and a right-handed component. 



propagating through the medium with momentum k^^ — (lu,K), in the presence of a gravitational potential and a 
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magnetic field B. In a homogeneous medium, the dispersion relation and wavefunction of the propagating modes are 
determined from the linear part of the effective field equation, which in momentum space takes the form 



i^-m- Y^eff + ^lY. ■ B)^j = . 
Se// is the self-energy of the neutrino, which in the situation that we are envisaging takes the form 



(A.l) 



(A.2) 



where 

= &mat + , (A.3) 

with and 6^ being given by the formulas quoted in Eqs. (|]^) and (|3.^). stands for the velocity four- vector 

of the background medium. We adopt the frame in which the background is at rest, so that = (1,0)- 

The formula for ha depends on the nature of the right-handed component field. If it is the conjugate field of a 
standard left-handed neutrino, then hj^ = — 6l. However, for a weak singlet field, which we assume to be the case in 



what follows, 6_R = 0. In Eq. (A.2) we are neglecting terms of order f/^/m^, as well as the magnetic-field dependent 



term that arises from the effective electromagnetic interactions of the neutrino in matter ||17|. The latter contribution 
is not significant for the values of the magnetic field that we are considering here. In addition, we have dropped the 
purely gravitational term denoted by hg in Rcf. which has a universal value for all the neutrinos (including the 
sterile ones) and therefore is not relevant for neutrino oscillations. 
Writing 

■ i 

V 



(A.4) 



in the Weyl representation of the gamma matrices, Eq. (A.l) becomes 



{uj-bi+ff- K)r] - TO^ + ■ B)^ = , 
(tj — (T • K)^ ~ mrj + fi{(7 ■ B)ri ~ , 



(A.5) 



where a are the Pauli matrices. Consider first the case in which B — 0. In this case, the equations have non-trivial 
solutions only if ^ and rj are proportional to the same spinor 0^ with definite helicity. The positive energy solutions 

2 2 

are then cj_ ~ + 6^ + ^ and cj+ ~ K + jj^ for A = — 1 and A = +1, respectively, with the corresponding Dirac 
wavef unctions 



where ?/, y' are of order m/K , and 



tp- = i>L + y 



ip+ = i^R + y' 



When i? 7^ 0, it is convenient to decompose B according to 

B = B^^K + Bt ■ 



(A.6) 



(A.7) 



(A., 



Then it is easy to see that, if Bt 7^ 0, the two helicity spinors get mixed in Eq. (A.5), and therefore we put 



?; = aq)^ + ecp 
e = P<f>+ + e'q 



(A.9) 



Substituting Eq. ( |A.9| ) into Eq. (A^) it follows that e and e' are of order m/K and fiB^/K, while the equations that 
are obtained for a and f3 can be recast in the form 



(A.IO) 
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where 



K + br 



H 



(m+M-B||)'' 
2K 



K + 



2K 



(A.ll) 



For a given solution of Eq. (A. 10), the Dirac wavefunction of the corresponding mode is obtained from Eqs. (A.4) 
and 



). Neglecting terms proportional to e and e', it is then given by 



(A.12) 



where ipL.R are defined in Eq. (A. 7). For an inhomogeneous medium, H is taken as the Hamiltonian for the amplitudes 

a, (3. 

The generalization to two or more families is straightforward. The same procedure yields an equation for the 
amplitudes that is identical to Eq. (A. 10), but where each element of H is considered to be a matrix in the neutrino 
flavor space. Similarly, the amplitudes a and /3 are spinors in that flavor space, and they determine the spin and 



flavor content of the corresponding propagating mode according to Eq. (A.12). These are the results quoted in Eqs. 
(ph and Q. 



In the second case that we consider, specified by the interaction Lagrangian terms given in Eqs. (3.12) and (3.13) 
the field equations, including the matter effects, are 



i)KR- + ■ B)UeL 

^KL-m^.^R 



0, 
0, 
0, 
0. 



where ill 



(A.13) 
(A. 14) 

Restricting ourselves for the moment to the situation in which _B|| =0, the positive energy solutions of these equations 
are found by putting, in the Weyl representation, 





"fl - °mat 



Veh 




NeR 


-{ 


Kr 




Kl 


- '-'( 



(A.15) 



Substituting these forms into Eq. (A.13), this procedure yields 



2K 

m 



2K 



—8- 



while the equations for a^^ and /3^^ can be written in the form 



with 



K 



H = 



"mat + °G 



2K 



-fiBj 



(A.16) 



(A.17) 



(A.18) 



mat 



When B\\ 7^ 0, the components of the wave function shown in Eq. (A.15) acquire an admixture of the opposite 
helicity spinor. Allowing this and carrying through the same steps as before, the resulting equation for a^^ and /3i7^ 
can be written once again as in Eq. ( |A.17| ), but with H given by the formula quoted in Eq. ( |3.15| ) . 
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